Optical response for the d-density wave model 
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We have calculated the optical conductivity and the Raman response for the d-density wave model, 
proposed as a possible explanation for the pseudogap seen in high T c cuprates. The total optical 
spectral weight remains approximately constant on opening of the pseudogap at fixed temperature. 
This occurs because there is a transfer of weight from the Drude peak to interband transitions 
across the pseudogap. The interband peak in the optical conductivity is prominent but becomes 
progressively reduced with increasing temperature, with impurity scattering, which distributes it 
over a larger energy range, and with inelastic scattering which can also shift its position, making it 
difficult to have a direct determination of the value of the pseudogap. Corresponding structure is 
seen in the optical scattering rate, but not necessarily at the same energies as in the conductivity. 



I. INTRODUCTION 

The pseudogap in the cuprates is widely believed to 
hold an important clue to the fundamental nature of 
these systemsii Through angular resolved photoemission 
studies 2 ^ of the variation around the Fermi surface of the 
leading edge of the electron spectral density, it has been 
established that the pseudogap has d-wave symmetry as 
does the superconducting gap. A smooth evolution of 
the superconducting gap into the pseudogap is also ob- 
served in tunneling 4 which shows that, in the underdoped 
regime, the gap fills in with increasing temperature but 
does not close as one goes through T c . These ideas have 
led to a school of thought that believes that the pseu- 
dogap is a precursor to superconductivity. One model 
is the preformed pair 5,6 model with phase coherence lost 
at the superconducting T Cl but the pairs remain till the 
higher pseudogap temperature T* . Another model in- 
cludes finite momentum pairsi, which exist at any finite 
temperature below T* but do not contribute to the super- 
fluid stiffness. No zero momentum Cooper pairs remain 
above T c where superconductivity is lost. 

A second school of thought envisages that the pseu- 
dogap has its origin in a competing order parameter^ 
The DDW theory is found among this group. It es- 
tablishes a charge density wave order which has d-wave 
symmetry and which opens at the antiferromagnetic Bril- 
louin zone as the competing order parameter. The DDW 
m0{ j e i8 4 9 4 io i ii 4 i2 breaks time reversal symmetry because 
it introduces bond currents which also have small or- 
bital magnetic moments. The two dimensional CuOi 
Brillouin zone is halved into a lower and an upper anti- 
ferromagnetic Brillouin zone. The DDW model can be 
included in the group of models associated with uncon- 
ventional density waves (UDW), which have been pro- 
posed to explain some phases of correlated electron sys- 
tems. For instance, UDWs have been used to under- 
stand quasi- two dimensional organic conductors, such as 
a-(BEDT-TTF) 2 MHg(SCN) 4 , with M = K, Rb and 
2^ ( i3 i i4 u nconven tional density waves have a gap which 
averages to zero on the Fermi surface (although in con- 
trast, the DDW gap, not opening at the Fermi surface, 



can only average to zero over the Brillouin zone). Like- 
wise, there is no periodic modulation in the charge or 
spin density, which makes them difficult to detect exper- 
imentally. For this reason it is very important to have 
predictions for transport properties in order to look for 
features that may characterize them. 

In this paper, we calculate the optical conductivity and 
the electronic Raman response for the DDW model. For 
tight-binding bands with first-nearest-neighbours only, 
the chemical potential acts as the lower cutoff on the 
interband signal in both optical responses as first dis- 
cussed by Yang and NayakiS in the case of optical con- 
ductivity. As the second-nearest-neighbour hopping is 
increased from zero, the cutoff remains, but shifts and 
is no longer at twice the chemical potential. The role 
of the chemical potential in the optical response is inti- 
mately related to the fact that the DDW order param- 
eter opens at the antiferromagnetic Brillouin zone. As 
a consequence optical spectral weight is transferred from 
intraband to interband transitions in the energy range 
of the gap with little change in total spectral weight. 
We obtain a peak in the real part of the conductivity 
coming from the interband contribution. This feature is 
robust with respect to different dispersion relations with 
second-nearest-neighbours. It is modified by impurities 
and by inelastic scattering, which shift the energy of the 
structure and obscure the cutoff, making an unambigu- 
ous determination of the pseudogap difficult. As the tem- 
perature is increased, the interband spectral weight loses 
intensity. It is fully depleted at the pseudogap tempera- 
ture T* at which point it has all been transferred to the 
intraband. A similar peak is found in the optical scat- 
tering rate, but it can be shifted in energy. This feature 
can be used to test for the DDW state, but so far, it has 
not been seen in the cuprates. 

The paper is organized as follows. In the next section, 
we present the DDW model, then in section ITTll we dis- 
cuss the one-particle properties relevant for our problem. 
Section lTvl treats the problem of the optical conductivity 
and of the optical spectral weight. Section El deals with 
the electronic Raman response. In section fvTl we con- 
sider the effects of elastic and inelastic scattering, and in 
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the final section our conclusions are presented. 



II. HAMILTONIAN 

We consider a two dimensional system of electrons with 
the Hamiltonian given by: 



H = 



1 



Ok - M)cL c k 



kk'<7<7' 

q 



^(k,k',q)4_ 



/Ck'cr'Cko-, (1) 



where c k and c kiCT are the creation and annihilation 
operators of electrons with momentum k over the entire 
Brillouin zone of the CuOi plane (2D) and spin a, fi 
is the chemical potential, e k is the dispersion relation 
and V(k, k',q) is the interaction between the electrons. 
Following Dora et alm^, the interaction is taken to be 
of the form F(k,k',q) = 2V/(k)/(k')<J(q - Q), where 
Q = (vr, it) is the wavevector where the density wave 
order takes place. 

Using the definition of the density wave order param- 
eter, 



where ± refers to the upper and the lower antiferromag- 
netic Brillouin zone, respectively. The imperfect nesting 
dispersion e k does not enter into the coherence factors 
in Eq. (JSJ, as has been previously discussed^, however, 
the effects induced by the imperfect nesting term enter 
in Eq. ©, where the it is outside the square-root. As a 
consequence, in all physical quantities, this term enters 
linearly along with the chemical potential. This fact will 
be of importance for the optical response. 

From the Hamiltonian of Eq. Q , we can write straight- 
forwardly the Green's function in Nambu notation, 
G(k, iujn)^ 1 = iu> n — Hjyjyw, a s follows: 

G(k, iw„) _1 = (iw n + (i - e+)f - e k f 3 + A k f 2 , (7) 

where the u> n are the fermion Matsubara frequencies. The 
spectral functions j4y-(k, u>) = — 2ImGij(k, u) + i5) are the 
following: 



Au(k,w) 

^22 (k, UJ) 

A 12 (k,uj) 



2tt{u15(uj + M - E+) + vIS(cj + 
2*{i£6(u + fi-E+) + u15(lu + /i - ££)}, 
A 3 i(k,w) 

2TTVkUk{S(uj + (i — E£) - 5(lu + fi — £7 k )}, 

(8) 



*A k = -VfQc) £ /(k')(cL +QCT Ckv), 

k 



(2) 



then, to within an additive constant, one arrives at the 
following effective Hamiltonian: 



H 



DDW 



E 



We observe that each Ay depends on the upper and lower 
bands, which reflects the physics of DDW theory where 
the quasiparticles spend time in both bands (i.e., the 
upper as well as the lower antiferromagnetic Brillouin 
zone) . 

We study the case with A k = ^(cos(k x ) — cos{k y )) and 

e k , where 



f(e k - n)cl a c ka + iA k 4. CT c k+Q(T . (3) the dispersion relation is of the form e k = e 



To diagonalize the Hamiltonian in Eq. we write it 
with momentum restricted to the first antiferromagnetic 
Brillouin zone using Nambu notation: 

Hddw = E *i( k )( e k^3 + (4 - V)r - A k f 2 W(k), 



(4) 

where ^(k) = (c^, a , c k+ q a ), $i stands for the Pauli 
matrices, e k = (e k — e k _q)/2 is the nesting dispersion 
relation and e k = (e k + e k _Q)/2 is the imperfect nesting 
one. 

Using a standard Bogoliubov transformation, the co- 
herence factors are found to be: 




e k = — 2t(cosfc r -|-cos/cj / ) and e k = — At' cos k x cos k y . We 
choose t' < 0, which is the relevant case for the cuprates. 
All the quantities are in units of the nearest-neighbour 
hopping parameter i, typically chosen to be between 0.1 — 
0.25eV. Here, t' is the next-nearest-neighbour hopping 
parameter. 



III. BANDS AND DENSITY OF STATES 

In the following we will comment on the band struc- 
ture given by Eq. JBJ together with the density of states. 
These one-particle properties can facilitate our under- 
standing of the optical response. In Fig. ^ we show the 
bands given by Eq. © and in Fig. [21 we plot the corre- 
sponding density of states given by: 



p(w) = [< 5 ^ - E £ ) + <S(oj - E k ) 



(9) 



wk 



(5) 



where -E k = \/ (e k ) 2 + A k . The energy eigenvalues are: 



Et = e+ ± Eu 



(6) 



where the momentum sum is over the entire Brillouin 
zone, and this quantity does not depend on the chem- 
ical potential as it can be absorbed into the frequency. 
First, consider the case where t' = 0. The upper, E£, 
and lower, EZ, bands are shown in Fig. 2Jl for A = 0.8t, 
illustrating that the gap is maximum at (tt, 0) and closes 
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at (7r/2,7r/2). For this band structure, at each value of 
momentum there are two branches, one for the upper and 
one for the lower antiferromagnetic Brillouin zone. This 
allows for vertical transitions with zero momentum trans- 
fer to be induced by the photon field from the lower to the 
upper band (interband transitions). Indicating schemat- 
ically a Fermi level /j on the band structure diagram in 
Fig. such a transition from the Fermi level at /_i in the 
lower band (hole doping) to the upper band, by an energy 
equal to 2\fi\, is shown. The significance of this result is 
that, in the absence of other mechanisms for absorption, 
such as impurity- or boson-assisted scattering, the mini- 
mum energy for an optical transition at finite frequency 
is 2|/x|. This 2\/i\ becomes the cutoff of the interband 
optical signal at low frequencies as will be seen from our 
full calculations. The corresponding density of states for 
this bandstructure is shown in Fig. [5] as the dashed line, 
which is the conventional one for a d-wave order param- 
eter. It has zero states at oj = 0. For this case of t' = 0, 
an analytical expression is available for the density of 
states in the continuum limit: = -K(^), where 

K(x 2 ) is the complete Elliptic integral of the first kind 
and N(0) is the density of states at the Fermi level in 
the state without the pseudogap. With t' nonzero, the 
bands and the density of states are modified. For small 
t' {\t'\ < A/4), the band structure is similar to the one at 
t' = with the lower band narrower and the upper band 
wider. If, on the other hand, t' is big (\t'\ > A/4), the 
bands overlap, as shown in Fig.^ 3 f° r t' = — 0.3< and the 
density of states is finite everywhere within the extended 
band. 



IV. OPTICAL CONDUCTIVITY 

We consider now the real part of the optical conduc- 
tivity <r(n). In the long wavelength limit (q — > 0) the 
current operator in momentum space is: 



TQ 



Ve k f 3 )* .(k, iuj n 



(10) 



where the £l n are the boson Matsubara frequencies. 
Within linear response theory, the optical conductivity 
is obtained from the current-current correlation function 
n(q -> 0,ifi„) as Re{cr QQ (f2)} = -(l/fi)Im{II a!Ce (i0 11 -> 
fl + iS)}, where a is a Cartesian coordinate, x, y, z. On 
the imaginary axis, Ii aa {iVt n ) is given in the bubble ap- 
proximation in terms of the Matsubara Green's functions 
G(k, iu> m ) by the following equation: 

II(ifi n ) = e 2 T ^ Tr[G(k,iw m +in„)(Ve^f + Ve k f 3 ) 



xG(k,iw m )(Ve k "f + Ve k f 3 ) 



(11) 



When the Green's function, given by Eq. Q, is expressed 
in terms of the spectral functions Ay (k, u>), the real part 
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FIG. 1: Energy dispersions E + and E of the DDW model 
for (a) t' = and (b) t' = -0.3/;, in the case of A = 0M. The 
horizontal lines indicate half filling (dotted) and a chemical 
potential 2\fi\ = l.Ot (dash-double-dotted, shown only in (a)). 
The tight-binding energy dispersion (dashed line) is presented 
for comparison. 




FIG. 2: Density of states for a DDW system with A = 0.8i, 
for t' = and t' = — 0.3t. The lower band narrows with in- 
creasing t and the density of states acquires significant finite 
weight within the extended band. 
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of the optical conductivity takes the form: 
duj f(w) - f(w + n) 



Re{a xx (r!)} = e 2 

£{[( 



2tt 



( <9/, , 



Aii(k,w + fi)Aii(k,a;) 



dk x > 

[(^) 2 -(l£)>-( k ^ +n )^( k ^)}' ^ 



where the Fermi function is given as f(u>) — l/(exp(/3w) + 
1) and (3 = l/fc^T. If we do not include disorder in the 
system, we can express the conductivity as Re{<r(f2)} = 
<To5(fi)+<Tjrit eT .(f2), where <to indicates the intraband con- 
ductivity or Drude weight and <Tmt er (n), the interband 
conductivity. Replacing the Aij by their expressions in 
, we obtain for the weight oq of the intraband conduc- 
tivity (Jiritra = ao5(Q): 



(To = ^ne* 



de+\ 2 e£/0ej 
£ k V dk. 



and for the interband conductivity for 17 > 0: 

- 2 sinh(/30/2) 

cosh/3(e+ - Ai) + cosh(/3f7/2) 



Winter (^) 



7re x - 



'der\ 2 A 2 



k ' - 2^ k ). 



(14) 



Concentrating first on the intraband conductivity given 
by Eq. I|13() . if we take both £' = and the continuum 
limit, we can extract analytical information. For t' = 0, 
Eq. ljT3*|) reduces to: 



h to 1, Eq. (|17|l can be finally stated in terms of Elliptic 
functions as follows: 



9i M 



2N(0)v F /w 



/ A 2 \ A 

£ (<^)- fori — ' 1 



{( 



{u + ^l) 2 ) 1 \w + h\ 
\u + n\ A ^ K f(Lu + f i) 2 



A 



+ 



A 2 



N + Ml 



for 7 



A 



/ (cj + /i) 2 \-| 
V A 2 // ' + 



> 1 
(18) 



where E{x 2 ) is the Elliptic integral of the second kind. 
In the limit of small and T relative to A, one can find 
that: 

(To ~ ^ , for /i = 0, and <To ~ f° r T = 0. 

(19) 

For /i = and T = 0, the Drude weight is zero, as 
expected from the density of states. As the temperature 
increases, spectral weight is transferred to the Drude, and 
likewise for finite \x. 

Next, we return to the case of the interband conductiv- 
ity. Setting t' = in (|14(l , the factor with the hyperbolic 
functions may be moved outside the momentum sum: 



Winter (^) 



sinh(/3ft/2) 



fl cosh(/3/i) +cosh(/30/2) 
where we have defined: 



y inter 

(ft), (20) 



k x k 

It is now possible to obtain an analytical expression for 
9inter(ty m the continuum limit: 



N(0) /2A\ 2 r /2 



cos" 



<l-(^) 2 cos 2 , 



CO 



^(o;) (15) which can also be written in terms of Elliptic functions: 



where 



Qinter 

(fi) 
t£jV(0)/7T 



5 1 M = E(|| 1 ) 2 J(^+' u -^)+^+M+^k)), 

k x k 

(16) 

which can be further reduced in the continuum limit, 
after integration over energy, to: 



M 



2 4H 



4A 2 



E 



4A 



* 2A 1 
for— < 1 



4A 2 



o9} 

ft 2 \ 1 . 2A 



E 



4A 5 



}' for 7T >1 - 

(23) 



ff iH = (^) 2 tv(o)- 



, ,4 r /2 



A 2 



(w + /i) : 



■ cos 2 9 



where we have made the angular dependence of Ag ex- 
plicit in Eq. i|17[l . Using the result that in two dimensions 
v Fx = v f/2i where vf is the Fermi velocity, and setting 



For small f2, QinteriP) — N(0)v F -g£ and for fl — > oo, it is 
equal to -^(0)t!f>(— ) 2 . Thus the interband contribution 
to the conductivity varies as (3VL at low frequency, with 
inclusion of the thermal factors, and as fl~ 3 at large Q, 
with a singularity expected at fl = 2 A. 

The function cosh( ff{f "Sn/2) ■ which multiplies 
<?inter(^) in Eq. (|20J) . is a universal function that shows 
the dependence of the interband conductivity on temper- 
ature and chemical potential. Taking the limit of zero 
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temperature, it becomes a step function 6>(f2 — 2\/j,\),i.e., 
zero, if f2 < 2\/i\ and one, otherwise. At finite tempera- 
ture, this step function becomes smeared. The fact that 
the interband conductivity starts abruptly at 2|/i is an 
important feature of the DDW modeli^. This feature 
does not arise, for example, in the case of the d-wave 
superconductor order parameter^, rather it is peculiar 
to the DDW model because the unconventional density 
wave gap opens up at the antiferromagnetic Brillouin 
zone boundary, while the d-wave superconductor order 
parameter opens up at the Fermi level. 

The interband conductivity for different values of the 
chemical potential is shown in Fig. [21 for a fixed value of 
the pseudogap, A = 0.8t, and temperature, T = 0.016i. 
Note that all figures for the conductivity are in arbitrary 
units. Fig. has been calculated using the continuum 
limit given by Eq. J5DJ and Eq. J33JI, and Fig. 03 using 
Eq. (|T4"|) with a lattice size 900 x 900. The main difference 
between the continuum limit and the discrete sum is that 
the singularity around 2A is reduced for the calculation 
performed by summing over the full Brillouin zone. Be- 
cause ai n ter(^) varies as /3f2 when O — > 0, the slope is set 
by 1/T, which becomes very steep when T is small, as 
seen in Fig. [3] (dotted curve). On the other hand when 
fi is finite and T — > 0, as we have already noted, the 
thermal factor in Eq. (|20[) provides a lower cutoff at 2|/x|, 
which is smeared in Fig. 3 due to finite temperature. 




Next, considering the effect of t' on the interband 
conductivity, we return to Eq. 1)140. The expression in 
Eq. (|14|l depends on t' only through the hyperbolic cosine 
where the imperfect nesting dispersion adds directly onto 
the chemical potential. Hence, the non-nesting term will 
alter the lower cutoff on the interband conductivity. In- 
deed, numerical calculation shows that the cutoff, previ- 
ously at 2|/i|, is shifted to lower frequency. The <7j n t eT .(fi) 
for different values of t' is shown in Fig. which is to 
be compared with frame (b) of Fig. [3J In all cases, we 
have taken half filling, maintaining the doping fixed. The 
doping x was defined, for temperature T, by: 



l-x = 2 



diop(uj)f(uj - fi), 



(24) 



where p(u>) is the density of states of Eq. JJjJ, shown in 
Fig. [21 which is normalized to one state over the entire 
Brillouin zone. The shape of the interband region is not 
much affected by the value of t' , but the cutoff at low en- 
ergy is now shifted down from 2|/x|, as discussed above. 
The surprising fact, that the shape of the interband con- 
tribution to the conductivity is not strongly affected by 
if , may be understood by referring to Fig.^ O n compar- 
ison of frames (a) and (b), the essential feature is that 
both bands shift at any value of k by the same amount , 
without significant change in the shape, which therefore 
does not affect the relative energy for vertical transitions. 



FIG. 3: Interband conductivity for t' = 0, A = 0.8t and 
T = 0.01 6t for (a) the continuum limit and (b) summing over 
all the Brillouin zone. The interband conductivity is shown for 
different values of the chemical potential jj,. In both figures, 
it can be seen that the interband conductivity rises abruptly 
at 2\fj,\, and at 2 A there is a singularity after which it drops 
rapidly. In the discrete case the singularity is reduced. 



We now turn to the total spectral weight, which has 
been discussed extensively as a key to understanding the 
mechanism of superconductivityiiiSiiS Having obtained 
expressions for inter- and intraband conductivities, it is 
straightforward to evaluate the different contributions to 
the total spectral weight, Wi n t ra = &o and Winter = 
J_ d£lo inter (^) ; & n d determine how they evolve with 
the various parameters of the model. 

Integrating expression i|14[) with respect to fi, the in- 
terband spectral weight is given by: 



Wi, 



E 

k 



dk a 



f{E--n)-f{E+-n) . (25) 



Thus, adding Eq. (|2*5^l to Eq. (fT3f (the intraband spec- 
tral weight) , we obtain for total spectral weight Wtotal = 
Winter + o"o . In Fig. [S] it is shown that the missing spec- 
tral weight, in the energy range from < 51 < 2|/i| in 
the interband conductivity, is transferred to the Drude 
weight as \fi\ increases. The solid curve applies to the 
case t' = 0, which corresponds to the dashed line in 
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FIG. 4: Interband conductivity for x — 0, A = 0.84, and 
T — 0.0164, for different values of t'. With t', the interband 
conductivity has a similar shape as for the case with t' = 0, 
but the lower cut off is shifted downwards in frequency. 




u7t 



FIG. 5: The interband and intraband spectral weights nor- 
malized to the total spectral weight for t' = and t' = — 0.3i, 
for A = 0.84. A transfer of spectral weight from interband to 
intraband occurs as the absolute value of the chemical poten- 
tial increases. 



Fig. [21 for the electronic density of states. For p = 0, 
all the optical spectral weight is in the interband contri- 
bution since p(u>) = at uj = 0, in this case, and there 
can be no intraband conductivity as a result. As \p\ is 
increased the intraband part increases rapidly and nearly 
linearly as does p(u) versus u of Fig. [3] in this range. As 
\p\ approaches A in value the growth in W intra becomes 
much less rapid and eventually saturates to being the 
total spectral weight. Thus, there is a transfer of spec- 
tral weight from the interband region to intraband when 
the doping is increased. While the results in Fig. [S] arc 
based on the lattice sum expressions for a 'intra an d 0"j nter , 
the results just described can be understood simply from 
Eq. (JTSJl in the continuum limit. At T — > 0, 
becomes a delta function and the integral in Eq. I|15fl be- 
comes equal to E(^ T ). This Elliptic function which de- 
scribes the dependence of the intraband spectral weight 
on p agrees remarkably well with the corresponding solid 
curve of Fig. [3] Next we consider the case for t' ^ 0. 
The dashed curves in Fig.0apply to the case 4' = —0.34, 
which corresponds to the solid curve in Fig. [3 where we 
see that there is now a finite density of states even at 
p = and hence the intraband contribution starts from 
a finite value of approximately 0.3 of the total weight. 
Because of the particular band structure involved, the 
intraband contribution first remains relatively flat until 
about p = —0.44, which is the value of the chemical po- 
tential at half filling for the particular parameters used 

here, A = 0.84 and 4' = -0.34. Beyond p 0.44, there 

is a linear behaviour similar to the one found for 4' = 0, 
although with a different slope, which finally saturates to 
where all the spectral weight is intraband. 

In Fig- El we show results for Winter (dashed), Wintra 
(dot-dashed) and W to tai (solid) as a function of the DDW 
gap A in units of 4. In these calculations, t' = and the 
doping was kept fixed at x = 0.1. Remarkably, the total 
spectral weight Wtotai is only decreased from its A = 



value by about 6% at A = 0.84. As A increases, the 
amount in the intraband component (which contains all 
the weight for A = 0) gradually decreases, however, it is 
still about twice the amount of that in the interband at 
the largest value of A considered. 
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FIG. 6: Total, interband and intraband spectral weights for 
fixed doping x = 0.1. The total spectral weight is not changed 
much by the opening of the DDW gap. 

To conclude, we wish to emphasize that our most im- 
portant results for the total spectral weight are: 1) there 
is a transfer in the spectral weight from the interband 
conductivity at low frequencies, < 51 < 2\p\, to the 
Drude weight at ft — 0, and 2) even though a pseudogap 
has opened, the total spectral weight is changed only by a 
few percent at fixed doping in the cases we have studied. 

Other aspects of the total spectral weight were con- 
sidered recently by Benfatto et al. 19 in relation to the 
temperature dependence of the total spectral weight and 
its relevance to the experimental results by Molegraaf 
et aliii, who found a very small decrease of Wtotai with 
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increasing temperature. Benfatto et al.'s work also ex- 
amined the effect of different choice of current operators. 
Aristov and ZeyherSfl have also briefly commented on the 
spectral weight issue with respect to temperature varia- 
tion. A feature of a recent work of Benfatto et aZiSi is 
the examination of the role of ward identities and vertex 
corrections on the temperature dependence of the optical 
sum rule, where they find that the total spectral weight 
increases with decreasing temperature. While there is 
some difference between their total spectral weight and 
ours, upon comparison with the frequency-dependent op- 
tical conductivity, we find very little difference between 
our calculation and the one shown in their paper. They 
are virtually the same, with only a small quantitative dif- 
ference at high frequency, and all the features discussed 
here remain intact. This indicates that these corrections, 
while important for the total spectral weight, may be ne- 
glected, as is common, for calculations of the frequency- 
dependent optical conductivity. Their work, done con- 
currently, is complementary to ours in that they examine 
primarily the issue of total spectral weight and we em- 
phasize the details of frequency-dependent conductivity. 
Our aim in this work is to determine how the total spec- 
tral weight is distributed in the optical conductivity and 
where it transfers to upon the opening of the pseudogap. 



V. ELECTRONIC RAMAN SCATTERING 

Now that we understand the crucial role of the chem- 
ical potential in the optical conductivity (indicating the 
frequency above which the interband response is different 
from zero), we examine its effects on the Raman response. 
The Raman response has the advantage that different re- 
gions of the Brillouin zone may be selected via the choices 
of incoming and outgoing photon polarizations. In con- 
trast, the optical conductivity is an average over all the 
Fermi surface. Hence, Raman experiments are relevant 
for unconventional density waves since they can provide 
information on the symmetry of the order parameter. 

In the point group D^h of the square lattice, we have 
three different channcls:£>i g , B 2g , and Ai g . The B\ g 
channel projects out the antinodal region of the Fermi 
surface, the B 2g , the nodal region, and the A\ g is a 
weighted average over the entire Brillouin zone. In eval- 
uating the Raman response for these different channels, 
we will discuss the continuum limit for simplicity. 

The real part of the optical conductivity multiplied by 
frequency, i.e. f2Re{er(S!)}, is related to the imaginary 
part of the Raman susceptibility^, x-y„ > m the case where 
the vp in the conductivity is replaced by j„, where 7^ is 
the Raman vertex for a particular channel. In the ab- 
sence of impurities, the Drude part of the Raman signal 
is zero and the interband part turns out to be: 



above: Ai g , B\ g and B 2gi and we have defined: 



Im{ X7 ^)} 



sinhQ3Q/2) 
cosh(/3/i) +cosh(/30/2) 
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xRc 



cos 2 6> 



(26) 



The expressions for the unrenormalized Raman vertices 
in the continuum limit are: 



7Ai e (0) = acos(46>), 
7Bi 9 (0) = &iCOs(2#), 
1B 2 M = b 2 sm(29), 



(27) 
(28) 
(29) 



and we have absorbed additional constant factors in the 

0.8r 




where v stands for the three different channels mentioned 
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FIG. 7: Raman signal for A\ g (dotted), B\ g (dashed), and 
B 2g (solid) channels for A = 0.5t, // = and T = 0.016*. At 
2A, there is a singularity in both A\ g and B\ g channels, and 
a small kink in B2 g . f2Re{cr(f2)} (dash-dotted) is also shown 
for comparison. All the signals are on the same scale. 



a and b's which are arbitrary. The integrals in Eq. I|26|l 
can be written in terms of Elliptic integrals which have 
the same form as for a d-wave superconductor^ The 
result is shown in Fig. for fi — and in Fig. 03 for 
/i = — 0.25<. To better illustrate the role of the chemi- 
cal potential in providing a lower frequency cutoff on the 
Raman signal, we have used a = b\ = b 2 = 1 • This is cer- 
tainly not realistic. In fact, in the tight-binding approxi- 
mation, the signal in B 2g is proportional to t' instead of 
and as a consequence we expect that the B 2g signal 
will be much smaller than the other ones. In the low 
frequency regime, the Ai g (dotted) and the B 2g (solid) 
signals are proportional to w/A, while B\ g (dashed) is 
proportional to (cj/A) 3 . For comparison we have also 
included our results for f2Re{cr(f2)} (dash-dotted curve) 
based on Fig. The Raman signal has recently been 
studied by Zeyher and Greco^ 5 . within the t — J model, 
with particular attention paid to the renormalization of 
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FIG. 8: Raman signal for the Ai g , Bi g , and B^g channels 
for A = 0.5*, fx = -0.25*, and T = 0.016*. The chemical 
potential cutoff is easily seen in B2 g and A\ g . The magnitude 
of the signal has been set to be the same in all cases, although 
Big will be, in reality, smaller. f2Re{<r(f2)} is also shown for 
comparison. 



the Raman vertices and to the role of superconductivity. 
Here, we are only interested in the mean field prediction 
for the region of the phase diagram where the DDW order 
prevails. 

From Fig. [SI we conclude that the effect of the 
cutoff is easily seen in the A\ g and the B^g signals, but 
less so for the B\ g channel, which is already quite small in 
magnitude at small lu. This is because the B\ g symmetry 
probes the antinodal quasiparticles where the pseudogap 
is largest. 



VI. SCATTERING 

We now turn to the problem of scattering by impuri- 
ties and inelastic scattering. We study only the optical 
conductivity since the results can be easily extrapolated 
to the Raman response. To include quasiparticle scat- 
tering, we replace the delta functions inside the spectral 
functions Ay(k, cj) in Eq. JSJ, with Lorentzians of the 
form: 



n(u + n-E±- Re£(k, lu)) 2 + 7 (w) 2 ' 

(30) 

where "/(lu) = — ImS(k, to) and E(k, uS) is the self-energy 
due to impurities or coupling to inelastic scattering. 
These Ay are inserted in the optical conductivity given 
by Eq. flfy . 

We study two cases, the scattering rate equal to a con- 
stant: y(u)) = 77, and the scattering rate of the marginal 
Fermi liquid: 7(0;) = 77 + 7rA|a>| where 77 and A are con- 
stants. 



A. 7(0;) = t) 

If 7(w) = rj, then ReS(k, lu) = 0. The effect of adding 
a constant scattering rate is that the interband contri- 
bution can no longer be separated from the Drude con- 
tribution. Still, these two contributions can be clearly 
discerned in the real part of the conductivity for a large 
range of parameters. For the numerical calculations of 




FIG. 9: Evolution of the real part of the optical conductivity 
with doping for a fixed gap, A = 0.8i, t' = 0, constant scat- 
tering rate 77 = 0.05*, and T = 0.016*. The values 2A and 2\fi\ 
are indicated by arrows. As the doping is increased, the lower 
edge corresponding to 2\n\ is moving to higher frequency and 
the Drude weight is increasing. 

the optical conductivity given by Eq. p2|l . we have used 
a lattice size of 300x300 for the Brillouin zone sum. The 
evolution with doping of the Re{o-(S!)} is shown in Fig. El 
In this figure, we can clearly see the redistribution of the 
spectral weight between the Drude contribution and the 
interband contribution when the doping is varied. In all 
these curves, solid for x — 0.04, dashed for x — 0.12 and 
dash-dotted for x=0.2, the lower edge of the interband 
region is at 2\n\ as indicated by arrows with some smear- 
ing due to temperature and impurity scattering. This 
edge moves upwards with increasing doping. There is 
little change however in the upper edge of the band at 
approximately 2A. For the dash-dotted curve, the in- 
terband region is almost completely gone and has been 
transfered to the Drude peak which has increased optical 
strength as compared to the solid curve. The effect of 
varying the impurity parameter 77 is very straighforward 
and, therefore, not shown. As r\ is increased the Drude 
component at low frequency is broadened in width and 
the peak drops, likewise the interband contribution is 
further smeared out and reduced in height. 

We now turn to the question of how Re{er(f2)} varies 
with temperature for fixed doping. For this calcula- 
tion we have used a fitted dependence of the gap with 
temperature^: A(T) = (1 - ±(T/T*) 4 )^1 - (T/T*) 4 , 
where we have chosen T* = 0.16*. As the temperature 
is increased, the gap seen as the upper edge of the in- 



9 



0.4 



a 

6 0.2 
"57 



o.o 



T=0.95T*,n=-0.262t 
T=0.85T*,n=-0.339t 
T=0.65T*,n=-0.419t 
T=0.1T*,n=-0.545t 




model: 26 ' 27 



1 

£2/t 



FIG. 10: Evolution of the real part of the optical conductiv- 
ity with temperature for A = 0.8i, i! = 0, T] = 0.05t and 
x — 0.08. The edge at 2 A is moving toward lower frequency 
with increasing temperature because of the temperature de- 
pendence of A. The Drude peak is increasing with tempera- 
ture, although the opposite behaviour is found when inelastic 
scattering is included. 



tcrband transitions moves to lower frequencies because 
the pseudogap closes. Secondly, the depression between 
low frequencies and 2|/i is filled in with increasing tem- 
perature. More importantly, the Drude weight increases 
with increasing temperature because of the collapse of the 
interband region. The DC conductivity displays a typi- 
cal semiconducting behaviour with resistivity decreasing 
with increasing temperature. However, in these calcula- 
tions we have kept the scattering rate constant. In any 
realistic case, the inelastic scattering will increase with 
T. This can be modelled by a temperature dependent 
scattering rate which increases with temperature, there- 
fore reducing the DC conductivity, and as a result there 
can be a change of its behaviour from semiconducting 
to metallic. Such inelastic scattering would modify the 
result of Fig. It would broaden the Drude peak rela- 
tively to what is seen and fill in the region between the 
intraband and interband contributions even more, ob- 
scuring the 2|/x| cutoff. 

Experimentally, the spectral weight distribution at low 
frequency and its temperature dependence has been stud- 
ied in reference^. For their underdoped sample, the 
spectral weight in the pseudogap region, does not move 
to lower frequencies with increasing temperature as is 
predicted in DDW theory, but rather it has the opposite 
behaviour; the spectral weight at low frequencies moves 
to the Drude with decreasing temperatureiiSiSii This be- 
haviour could only be obtained in the DDW model with 
strong inelastic scattering. 

Next, we calculate another important quantity that 
is obtained in optical experiments, namely the optical 
scattering rate. For this, we use the extended Drude 



n 2 



Rc{er(ft)} 



4tt Re{a{Q)} 2 +lm{a{Q)} 2 



(31) 



where we have calculated numerically the Im{a(£l)} by 
Kramers-Kronig relations, and for the plasma frequency 
we have used that derived from the Drude weight in the 
case of A = 0. The result is shown in Fig. ^] The most 
striking feature in this figure is a peak which arises from 
the interband contribution. This is an important result 
since, in principle, it can be verified experimentally. We 
note that as the temperature increases, the size of the 
peak decreases in magnitude and shifts to lower energy 
before disappearing completely at the pseudogap temper- 
ature T* = 0.16t. Also the position of the peak in the 
scattering rate of Eq. (|31[) is not at the same position 
as the peak in the real part of the conductivity. It is 
lower. In particular for the solid curve in Fig. 1111 which 
is the highest temperature shown, the peak falls below 
0.5i while the corresponding structure in the real part 
of the conductivity in Fig. ^| is nearer 0.8i, although in 
that case it cannot be characterized simply as a peak. 
Later, we will comment on the origin of these shifts, be- 
tween the scattering rate and the optical conductivity, of 
the interband structure. A second feature to be noted in 
Fig. ^] is that the scattering rates at u> — > are consid- 
erably larger than twice the value of the impurity scat- 
tering rate that we have used. This feature will also be 
addressed at the end of this section. Finally, as discussed 
before, if the impurity parameter r\ is increased, the con- 
ductivity is broadened and the effect on the peak seen in 
the scattering rate is for it to be broadened and reduced. 
Recently a peak in r op 1 (cj) in the region of 400 to 3500 
cm -1 has been observed in the cobaltates^. This peak 
reduces and eventually vanishes with increasing temper- 
ature. There appears to be a corresponding peak in the 
conductivity, as well, at higher frequency (400 to 6000 
cm -1 ), not unlike what we find in our model. 



B. Inelastic scattering 

To present a more realistic comparison with experi- 
ments, we should include the inelastic scattering. It is 
important to check whether or not the strong feature 
found above, in the optical scattering rate, will be washed 
out by the inelastic scattering. We can do this using the 
self-energy of the marginal Fermi liquid 29 in our spec- 
tral functions Ay(k, lu). This phenomenological model 
assumes that the excited charge carriers interact with a 
wide, flat spectrum of excitations over T < lo < lo c , where 
lu c is a high energy cutoff of the spectrum. This model 
was introduced to explain the observed linear behaviour 
in the scattering rate T ~ max(|a;|,T) measured in the 
optical conductivity, as well as in the resistivity, and in 
the Raman cross-section^ 7 ^ " The marginal Fermi liquid 
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FIG. 11: Evolution of the optical scattering rate with tem- 
perature for x = 0.08, A = 0.8*, t' = 0, and rj = 0.05*. 
At low temperatures, there is a strong peak caused by the 
interband processes, which washes out with increasing tem- 
perature, and the scattering rate eventually recovers the pure 
Drude behaviour. 
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FIG. 12: Comparison of the real part of the optical conduc- 
tivity including the marginal Fermi liquid self-energy (dash- 
dotted curve) with the case where the scattering rate is con- 
stant (dotted curve): j(u>) = rj. Here, A = 0.8t, fi = — 0.5i, 
t' = 0, ttA = 0.5, 77 = O.lt, and T = 0.016*. Two characteris- 
tics of the marginal Fermi liquid are evident: 1. there is more 
spectral weight at higher frequencies and 2. the low frequency 
peak is narrowed in comparison with the dotted curve. Also 
shown are the results with a DDW included. 



ansatz is: 27,29 



^mfl(w) = 2Xuj In wrAM, 



(32) 



where A is a dimensionless coupling constant. 

Results based on the MFL self-energy are shown in 
Fig. 1121 where they are also compared with previous 
results for impurities alone (i.e., a constant scattering 
rate). The dotted line is for r/ — O.lt and the dash-dotted 
curve is for 77 = 0.1* plus the MFL contribution to the 
self-energy. On comparing these two curves, both calcu- 
lated at temperature T = 0.016*, we note that consider- 




FIG. 13: Optical scattering rate for the MFL with the DDW. 
Shown is l/r(f2) (solid line) for the case of A = 0.8*, /1 = 
—0.5*, *' = 0, 7rA = 0.5 and r\ — 0.1*, calculated from the 
conductivity curve given by the solid line in Fig. 1121 The 
dashed line is an approximate fit to the conductivity curve of 
Fig. I12l using the model of a MFL plus a displaced oscillator, 
as discussed in the text, and from this 1/t*(Q) is calculated 
(dash-dotted line) and can be seen to be related to l/r(f2) of 
Eq. 131H by a multiplicative constant. 



able optical spectral weight is lost at small fl in the dash- 
dotted curve, as compared to the dotted Drude curve. 
There is a crossing at Cl/t slightly above 1 and for the 
MFL there is a long tail in the mid-infrared region of the 
optical spectrum, extending to high Q and staying well 
above the Drude curve. These tails are an essential fea- 
ture of the conductivity in the cuprates, an indication of 
a strong electron-boson inelastic contribution (Holstein 
band). The depletion of the optical weight at small 
reflects the fact that the quasiparticle spectral weight in 
the MFL, Z(u) = (1 - varies as (2Aln|^)- 1 

at small u> and therefore vanishes logarithmically at the 
Fermi surface. The other two curves in Fig. 1121 include 
a DDW. The dashed curve is based on a constant quasi- 
particle scattering rate of 77 = 0.1 while the solid curve is 
for 0.1 + 0.5 |w|. The interband contribution manifests it- 
self as a broad peak centered around fz ~ 1.4*, in the 
dashed curve, which is also reasonably well separated 
from the Drude contribution centered around O = 0. 
As the impurity scattering is increased, the region be- 
tween the interband and the intraband fills in, and the 
interband peak broadens and is reduced in height (not 
shown). In a sense, the MFL result (black curve) can 
be understood as a case where impurity scattering effec- 
tively increases with increasing frequency. This has the 
outcome of spreading the interband spectral weight over 
a much larger frequency range reducing the prominence 
of the peak seen in the dashed curve and shifting it to 
lower energy. Because of this shift, it is not possible to 
accurately determine the size of the pseudogap from the 
effective position of the interband transitions in the real 
part of the conductivity. Nevertheless, a peak still re- 
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mains in the MFL case, although now it is not so well 
separated from the intraband contribution at low ener- 
gies and the 2|/x| cutoff is no longer visible. 

We also note from Eq. 112|l . that for the DC conduc- 
tivity at zero temperature, the electron spectral func- 
tions Aij (k, lj) are all evaluated at zero frequency. In 
this case, only the impurity self-energy remains, since 
Y,mfl(ui = 0) = in Eq. 132|) . therefore, both the im- 
purity and the marginal Fermi liquid models will have 
the same DC value because we have chosen the same r\ 
and plasma frequency in both cases. However, for finite 
temperature, the marginal Fermi liquid DC conductivity 
will be reduced because of the inelastic scattering, a fea- 
ture not present in the impurity case. When the DDW 
gap is included, spectral weight is transferred to the in- 
terband processes. In both cases of impurity only and 
impurity plus MFL, this transfer of spectral weight is 
about 50% for the parameters chosen, namely /i = — 0.5£ 
and A = 0.8i. This is consistent with the solid curves in 
Fig. [S] which cross at this value of the chemical potential 
for t' — 0. For the impurity case the value of the DC 
conductivity is reduced to 53% of its original value while 
for the marginal Fermi liquid the percentage is 47%. 

In Fig. El we show our results for the optical scatter- 
ing rate obtained from the marginal Fermi conductivity 
of Fig. El What is shown as the solid curve is the usual 
1/t op (lu) of Eq. H31fl . It is this quantity that is used in 
the experimental literature. However, it is useful to de- 
fine a new scattering rate denoted by l/r* p (w) which is 
related to l/r op (w) by a simple multiplicative constant. 
For 1/ t* p ( u), what is used for forming the righthand side 
of Eq. Q31[) is the value of the DC conductivity divided by 
2rj, rather than flp/An. For the impurity case without the 
DDW, this would give Sip/An, as would also be the case 
for the pure MFL at T = 0. At T ^ 0, an inelastic con- 
tribution to the scattering would also need to be added 
to 2-q. For the case shown here, where T = 0.016t, this 
is about a 10% correction. As we have seen, when the 
DDW is present, considerable spectral weight is trans- 
ferred from the intraband to the interband absorption, 
and so the Drude centered about SI = is strongly de- 
pleted, and the DC conductivity is correspondingly re- 
duced. This means that l/r* p (w) is about one half the 
value of 1/t op {uj) for the case considered here when the 
DDW is present. An important feature of 1/t* p (uj) is that 
as u! — > 0, it agrees exactly with the residual scattering 
rate 2r] with two modifications. First, when quasiparti- 
cle scattering is present, the interband contribution can 
make a contribution to the DC conductivity, and second, 
at finite temperature in the MFL, there can also be a 
small increase in the DC scattering due to the inelastic 
scattering. Both these effects can modify somewhat the 
u) — > of the 1/t* (ui). Returning now to Fig. ^2 we 
can understand why, in the lower temperature run, the 
DC limit of 1/t op (u>) is more than twice the value of 2r\. 
It reduces with increasing temperature but even for the 
solid curve at T = 0.95T*, it is still above 2r\ = O.Lt. Re- 
turning to the solid curve of Fig. El we note that its DC 



value is also above 2-q mainly due to the depletion of the 
intraband contribution. At higher frequencies, it shows 
the interband peak slightly below f2/f = 0.5 and then 
takes on a quasilincar behaviour as expected of the pure 
MFL, but with a slope which is considerably reduced over 
the input scattering rate, shown as the dotted curve. The 
reduction in slope can be traced to the presence of the 
interband transitions as will be made clear in the follow- 
ing. 

We now turn to the shift in the interband structure to- 
ward lower frequency seen in the scattering rate of Fig. 1131 
when compared to the structure in the conductivity of 
Fig. El This also holds for Fig. ^2 when compared to 
Fig. El Some understanding of this can be obtained 
from the study of a simplified model. For convenience, 
we use a simplified form for the marginal Fermi liquid 
conductivity^ 



MFL 



MFL\2 



4tt 



-iu)(l ~ 2Alog \uj/2uj c \) + (277 + 7rAw)' 

(33) 

Note that for the conductivity, it is 2i] + ir\\u\ which 
enters as the damping rate. The constant impurity scat- 
tering rate is twice its quasiparticle counterpart, but the 
MFL piece remains unaltered. Eq. (|33|l is a reasonably 
useful approximation to the more accurate calculation 
that would arise from implementing the exact form for 
the normal state conductivity^ 



4tt v J Q v — S(o;) — Yi(y — u) + i2r/ 



(34) 



with the marginal Fermi liquid self-energy given in 
Eq. (|32[) . Adding to Eq. I)33|l a form based on a displaced 
Lorentzian oscillator in the dielectric function of energy 
loe, which is used to model the interband contribution: 32 



(O 



4-7T (to 



(35) 



we form the optical scattering rate defined in Eq. Q31(l 
for the combined system with one change. Instead of 
using the complete plasma frequency coming from inter- 
band and intraband combined, we use only the intra- 
band contribution. This is done to define 1/t* p (lo). We 
use Eq. I|33l) and Eq. (|35|) to model as well as possible 
the solid curve (MFL+DDW) in Fig. EJ The result is 
shown in Fig. El as the dashed line. The parameters 
used are A = 1/(2tt), 77 = O.lt, V = 0.9i, lo e = 0.6£ 
and {Q^or /n™ FL ) 2 = 0.57. With these parameters we 
can calculate the scattering rate l/r*(il) for our model. 
It is given by the dash-dotted curve which can be com- 
pared with the dotted straight line equal to 0.2 + 0.5fl/t. 
These quantities agree at f2 = as expected. It is to be 
noted that this graph provides a good qualitative confir- 
mation of the complete numerical results of Fig. 1131 fsolid 
curve). We note that l/r*(fi) differs from \/t(Q) of for- 
mula Eq. (|31|) mainly by a constant numerical factor. We 
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also note that the peak in the scattering rate has shifted 
downward as compared with its position in the conduc- 
tivity. This is clearly due to the particular combination 
of real and imaginary part of the conductivity defined in 
Eq. (|3 ll> . The reduced slope of the quasi- linear behaviour 
of the dash-dotted curve at large f2 when compared with 
the pure MFL slope is reproduced in our model (solid 
curve). It can now be seen that it has its origin in the fact 
that an interband contribution is present and is not due 
to an alteration of the MFL representation of the remain- 
ing intraband contribution (with the lorentzian term set 
to zero, we recover the straight dotted line as expected). 
We can conclude from this analysis that, to good approxi- 
mation, the complete conductivity can be nicely modelled 
as a superposition of a MFL component of reduced spec- 
tral weight, but otherwise unaltered, plus an interband 
component given approximately by Eq. (|35|l . This repro- 
duces the full numerical calculation reasonably well, how- 
ever, the original pseudogap value and cutoff 2|/x| could 
never be inferred from fitting the experimental data to 
this model. 

In principle, the formation of a DDW will also lead to 
changes in the self-energy of the charge carriers. For ex- 
ample in the DDW model the electronic density of states 
acquires additional energy dependence as seen in Fig. [21 
which, in turn alters the quasiparticle scatteringii&M 
Here, we have modelled this scattering either with a con- 
stant elastic scattering rate or, to be more realistic, by a 
marginal Fermi liquid (MFL) contribution. This contri- 
bution to the inelastic scattering is directly responsible 
for the quasilinear dependence, in frequency to, of the 
optical scattering rate, seen in Fig. ^| at u> beyond the 
interband contribution. This quasilinear dependence is a 
hallmark of the cuprates and is central to the MFL phe- 
nomenology which was introduced to correlate a large 
number of the normal state properties of the high T c ox- 
ides, which are anomalous. Unfortunately, there is no 
generally accepted underlying microscopic framework 35 
to the marginal Fermi liquid ansatz, which could be em- 
ployed to introduced the necessary modifications to the 
self-energy that result from the DDW formation. Such 
considerations go beyond the scope of the present study. 

The effects of energy dependence in the electronic den- 
sity of states p(e) on the self-energy and on the conductiv- 
ity have recently been studied 3 ^ for elastic impurity scat- 
tering and are well understood. For example, a depres- 
sion in p(e) below its constant background value around 
the Fermi energy, leads directly to a corresponding reduc- 
tion in the optical scattering rate in this same frequency 
region. In a non-self-consistent approach, the constant 
quasiparticle scattering rate of the flat p(e) case, is simply 
multiplied by p(e) to lowest order. In our case, this would 
be the p(e) of the DDW as shown in Fig. [21 This would 
correspondingly reduce the scattering rate of Fig. ED at 
small u>. However, it was found in reference 3 ^ that such 
a procedure overestimates the effect of p(e) on the quasi- 
particle as well as on the optical scattering rate. In a 
self-consistent approach, the impurities and, even more 



importantly, the inelastic scattering of the MFL smear 
out the structures in p{e) which are due to the opening 
of the DDW gap. Nevertheless, a reduction of r _1 (J7) 
below the value seen in Fig ll^l is expected for frequencies 
below the gap and would be associated with self-energy 
corrections. 

Finally, we discuss the issue of comparison with exper- 
iment. This model predicts the appearance of an inter- 
band contribution in the pseudogap phase and an accom- 
panying peak in the optical scattering rate. This inter- 
band contribution should have a feature at lower energy 
corresponding to a lower cutoff set by 2|/i| and a higher 
cutoff of about 2 A. The lower feature should shift to 
higher energy with increased doping, however, from ex- 
periment the gap decreases with doping and hence the 
upper edge in the conductivity should decrease with in- 
creased doping. This has not been seen in experiments 
on the high T c cuprates. Examples of optical experiments 
examining the underdoped state of various cuprates are 
given in R e f s , 26 i 37 i 38 for YBCO, Bi2212, and LSCO mate- 
rials. There is no evidence for an interband feature with 
energy scales as described here and the optical scattering 
rate2& does not show a peak. The most up-to-date review 
on optical conductivity experiments in high T c cuprates 39 
does not give further evidence that could support this 
model. Indeed the DDW model has been controversial 
in its applicability to experiments in the cuprates. Some 
criticisms have involved the variance of the model with 
regard to photoemission— and tunneling^i, but some 
of these issues have been possibly rectified by further 
calculation! 5 -^. Here we find that the experimental ev- 
idence from optical conductivity, at this time, does not 
support the DDW model for the cuprates. However, this 
model may find application in other systems in the fu- 
ture. 



VII. CONCLUSIONS 

We have calculated the optical conductivity and the 
Raman response for a DDW system. There is a cutoff 
at low frequencies in the interband response for both the 
optical conductivity and the Raman response. This cut- 
off is at 2 1 /it | if there is no imperfect nesting term in the 
dispersion relation. The imperfect nesting adds directly 
to the chemical potential and shifts the cutoff to lower 
frequency but does not eliminate it. Including scattering 
makes the cutoff harder to discern. 

As the chemical potential is changed we find a read- 
justment of the spectral weight between the inter- and 
intraband conductivities. That part of the interband con- 
tribution, ranging from frequencies zero to 2\/j,\, is trans- 
ferred at zero temperature to the Drude weight as the 
chemical potential is increased with increasing doping. 
However, contrary to a first expectation there is little 
loss in the total spectral weight with the opening of the 
pseudogap. This result agrees with experiments^ 

As mentioned above, the 2\fi\ cutoff also appears in 
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the electronic Raman scattering. The best channels to 
observe the cutoff at 2 1 /x | are B 2g and A\ g . The B lg 
channel projects out the antinodal region where the gap 
is biggest, and the low intensity of the Raman scattering 
at small u>, where it goes like (w/A) 3 , may make the 
cutoff difficult to see. 

The behaviour with temperature of the optical con- 
ductivity has two characteristics: 1) the expected shift 
of the 2A feature to lower frequencies, for a mean field 
theory with order parameter vanishing at T* , and 2) the 
transfer of spectral weight from interband to intraband 
as the interband contribution collapses. There is also a 
filling of the depression between low frequencies and 2\fi\. 

An important result of this paper is that there is a peak 
in the real part of the optical conductivity in the region 
Q c < fi < 2A that comes from the interband conduc- 
tivity with £l c = 2|/i| for t' — 0. This peak is robust 
against different dispersion relations including second- 
nearest-neighbours. It is reduced by the impurity scatter- 
ing and by the inelastic scattering, which can also shift 



it. It is washed out with increasing temperature van- 
ishing completely at the pseudogap critical temperature. 
Similar structure is found in the optical scattering rate, 
but it is shifted relative to that in the optical conduc- 
tivity. With strong inelastic scattering, as may be found 
in the MFL model, the 2|/i| cutoff and 2 A features be- 
come more difficult to determine directly. To the best of 
our knowledge, this interband structure has not yet been 
identified in experiments on cuprates. 
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